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Abstract 

We consider Feynman-Dyson’s proof of Maxwell’s equations using the Jacobi identities 
on the velocity phase space. In this paper we generalize the Feynman-Dyson’s scheme by 
incorporating the non-commutativity between various spatial coordinates along with the ve¬ 
locity coordinates. This allows us to study a generalized class of Hamiltonian systems. We 
explore various dynamical flows associated to the Souriau form associated to this generalized 
Feynman-Dyson’s scheme. Moreover, using the Souriau form we show that these new classes 
of generalized systems are volume preserving mechanical systems. 

Keywords : Feynman problem, Souriau form, noncommutativity, generalized Hamiltonian dynam¬ 
ics. 

PACS numbers: ll.10.Nx, 02.40.Yy,45.20.Jj. 

1 Introduction 

The study of exotic particle models with non-commutative position coordinates was started in the 
last decade. There are several physical phenomena appearing in condensed matter physics, namely 
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semiclassical Bloch electron phenomena, fractional quantum Hall effect, double special relativity 
models, etc., that exhibit such feature. All these models share the somewhat unusual feature that 
the Poisson brackets of the planar coordinates do not vanish. This class of dynamical structures has 
appeared in geometric mechanics and geometric control theory [Q1l2l 3] [4j . In her thesis, Sanchez 
de Alvarez fl5l indicates a characterization of the Poisson structure in terms of Poisson brackets of 
particular functions on the tangent bundle TP of a Poisson manifold P, and discusses its functorial 
properties. 

A very noble derivation of a pair of Maxwell equations was originally proposed by Feynman, 
but the exact details of his argument came to the scientific community from the work of Freeman 
Dyson According to Dyson, Feynman showed him the construction and examples of the 
Lorentz force law and the homogeneous Maxwell equations in 1948. A derivation of a pair of 
Maxwell equations and the Lorentz force is based on the commutation relations between positions 
and velocities for a single non-relativistic particle. In general the locality property that different 
coordinates commute is assumed. Due to increasing interest in non-commutative field theories, 
it is worthwhile to consider the non-commutative analogue of Feynman approach. This destroys 
the axiom of locality, which, according to Dyson, was the original aim of Feynman. Tanimura 
ITT] gave both a special relativistic and a general relativistic versions of Feynman’s derivation. 
Land et al. HO examined Tanimura’s derivation in the framework of the proper time method in 
relativistic mechanics and showed that Tanimura’s result then corresponds to the five-dimensional 
electromagnetic theory previously derived from a Stueckelberg-type quantum theory in which one 
gauges the invariant parameter in the proper time method. An extension of Tanimura’s method 
has been achieved [[9] by using the Hodge duality to derive the two groups of Maxwell’s equations 
with a magnetic monopole in flat and in curved spaces. A rigorous mathematical description of 
Feynman’s derivation connected to the inverse problem for Poisson geometry has been formulated 
in lUOl (see also 02). Hughes 02 considered Feynman’s derivation in the framework of the 
Helmholtz inverse problem for the calculus of variations (see also P3l and 04l f. 

In fact, it was pointed out by Jackiw that Heisenberg suggested in a letter to Peierls that spatial 
coordinates may not commute, Peierls communicated the same idea to Pauli, who informed it to 
Openheimer; eventually the idea arrived to Snyder [131 [16] who wrote the first paper on the subject. 
Nowadays the physics in non-commutative planes is relevant not only in string theory but also in 
condensed matters physics ifTTl . In the context of the Feynman’s derivation of electrodynamics, it 
has been shown that non-commutativity allows other particle dynamics than the standard formalism 
of electrodynamics |[T8l . Noncommutative quantum mechanics is recently the subject of a wide 
range of works from particle physics to condensed matter physics. This has also been studied from 
the point of view of Feynman’s formalism in ||T9l . 

The examples of exotic mechanics started to appear around 1995. Physicists obtained various 
models such that the Poisson brackets of the planar coordinates do not vanish. Souriau’s orbit 
method [201Efl 22] was used to construct a classical mechanics associated with Levy-Leblond’s 
exotic Galilean symmetry. In terms of the Souriau 2-form a wide set of Hamiltonian dynamical 
systems have been described in [24l [25 , .261 [27]. Levy-Leblond lf28l has realized that due to the 
commutativity of the rotation group 0(2, M), the Lie algebra of the Galilei group in the plane 
admits a second exotic extension defined by 

[K u K 2 ]=mI, 
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where k is the new extension parameter. For a free particle the usual equations of motions are 
unchanged and k only contributes to the conserved quantities. It yields the non-commutativity of 
the position coordinates. 

Feynman procedure to obtain Maxwell’s equation in electrodynamics has been reviewed under 
different kind of settings, and several nontrivial and interesting generalizations are possible, see for 
instance (29l[30l[3Tl[32l[33l[34l[35l[36l|. Recently, Duval and Horvathy |[29Tl successfully applied 
the techniques of Souriau’s orbit method lU71 l38,39l to various models. Incidentally, one of these 
models can be viewed as the non-relativistic counterpart of the relativistic anyon considered before 
by Jackiw and Nair ll32l . Mathematically, the ‘exotic’ model arises due to the particular properties 
of the plane. A wide set of dynamical systems can be derived from the Lagrange-Souriau 2-form 
approach in three dimensions and the generalizations to higher number of degrees of freedom have 
been outlined in f[26ll . 

Wong’s equations describe the interaction between the Yang-Mills field and an isotopic-spin 
carrying particle in the classical limit. Feynman-Dyson’s proof offers a way to check the consis¬ 
tency of these equations fiOil . See also iBTll for a very recent paper. 

In a slightly different context Kauffman ll42ll introduced discrete physics based on a non- 
commutative calculus of finite differences. This gives a context for the Feynman-Dyson derivation 
of non-commutative electromagnetism. More recently, Kauffman Il43ll found an interesting way to 
describe mechanics in a curved background interacting with gauge fields in such a way that the 
physical equations of motion emerge automatically from underlying algebraic relations in a non- 
commutative geometry and this construction depends largely on the Feynman-Dyson construction. 
In an interesting paper, Cortese and Garcia ll44ll studied a variational principle for noncommutative 
dynamical systems in the configuration space. In particular they showed that the non-commutative 
consistency conditions (NCCC), that come from the analysis of the dynamical compatibility, are 
not the Helmholtz conditions of the generalized inverse problem of the calculus of variations. 
It has been shown that the 0-deformed Helmholtz conditions are connected to a third-order time 
derivative system of differential equations. Noncommutative phase spaces have been introduced by 
minimal couplings in ||45Tl and then some of them are realized as coadjoint orbits of the anisotropic 
Newton-Hooke groups in two- and three-dimensional spaces. This has been further generalized to 
realize noncommutative phase spaces as coadjoint orbit extensions of the Aristotle group in a two 
dimensional space f46] . 

In this article we apply Souriau’s orbit method to study exotic mechanics on the tangent bundle 
or velocity phase space. Souriau first unified both the symplectic structure and the Hamiltonian into 
a single two-form. It has an exotic symplectic form and a free Hamiltonian and yields a generalized 
Hamiltonian mechanics. Duval and Horvathy used Souriau’s orbit method to construct a classical 
planar system associated with Levy-Leblond’s two-fold extended Galilean symmetry. The four 
dimensional phase space is endowed with the following exotic form 

Q 

Q = dpi A dqi + —€ij dpi A dpj, 

where summation on repeated indices is understood. The exotic term in the symplectic form only 
exists in the plane. Following [47,33 we also explore a volume-preserving flow on a symplectic 
manifold from the Souriau form associated with the velocity phase space. 

Many authors [29], [30, 49] have generalized this modification of the symplectic form by intro- 
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ducing the so-called dual magnetic field such that 


VL = dpi A dqi + -gij dpi A dpj + -/„ dq { A c%. 

The coefficients and f VJ are responsible of the noncommutativity of momenta and positions, 
respectively. The classical dynamics in noncommutative space leads to noncommutative Newton’s 
second law [ 50lf5T|| . This generalization can be studied in various types of noncommutative space- 
times; for instance Harikumar and Kapur studied in If52l the modification to Newton’s second law 
due to the kappa-deformation. In a very recent paper the modification of integrable models in the 
kappa-deformed scheme is analyzed [T53T1 and kappa-Minkowski space-time through exotic oscil¬ 
lator is studied in [f54l . Zhang et al. 155 . [561 studied the 3D mechanics with non-commutativity, 
where the potential may also depend on the momentum. They obtained the conserved quantities 
by using van Holten’s covariant framework. It is known that the Snyder model has the remarkable 
property of leaving the Lorentz invariance intact. Recently, motivated from loop quantum gravity 
an idea has been proposed to extend the Snyder model 031 to space-times of constant curvature, 
by introducing a new fundamental constant whose inverse is proportional to the inverse of the cos¬ 
mological constant. More recently, classical dynamics on Snyder space-times has drawn a lot of 
attention to physicists It57l[58lf59u 60[|. Moreover Snyder dynamics in curved space-time has been 
extended by Mignemi et al. in [f6TL l62ll . See also lf63l . 

The main theme of our paper is to show that non-commutativity between coordinates allows us 
to construct various other generalized classes of dynamical systems. Tools of non-commutative ge¬ 
ometry often appear in quantum gravity. Using a differential geometric theory on non-commutative 
space-time Aschieri et al. If64l defined ^-deformed Einstein-Hilbert action, and by means of their 
technique of the deformation of the algebra of diffeomorphisms one can derive ^-deformed inte¬ 
grable systems f[65l and Newtonian mechanics |[66l . Today we find non-commutativity in various 
fields of modern physics such as, graphene. Hydrogen atom spectrum, etc. Il67l [68,[69|. 

This paper is organized as follows: in order to the paper be self-contained, Section 2 is devoted 
to a review of multivector fields, Poisson bivector, Schouten-Nijenhuis bracket and various other 
geometrical tools. We give a brief geometrical description of Poisson manifolds in Section 3. Sec¬ 
tion 4 is devoted to Souriau’s formalism of generalized symplectic forms. We illustrate Souriau’s 
construction through examples. Section 5 is focused on the construction of Feynman-Dyson’s 
scheme and its connection to Souriau’s method. Section 6 relates volume preserving mechanical 
systems and Souriau’s form. We finish our paper with an outlook in Section 7. 

2 Geometrical background 

Let T(M) be the algebra of U^-class functions (the algebra of classical observables) on a manifold 
M (the classical state space). We denote by IP (A/) the space of C 00 -class differentiable p-forms, 
and by A P (M ) the space of 6 ,oc -elass skew-symmetric contravariant tensor fields of order p, often 
called p-vectors. By convention we set A°(M) = f2°(Af) = D(M) and A P (M) = fP(M) = 0, 
whenp < 0. Then, 

Q(M) = 0 fP(M) and A(M) = 0 A P (M), 

p&i 
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are Z-graded algebras under their exterior products; moreover both are anticommutative, so, for 
instance, if P £ A P (M) and Q £ A q (M) 

PAQ = -(-I)mqaP. 

When a is a 1-form and X is a vector field the C°°(M)-c lass function (a, X) given by 

(a,X)(x) := (a(x),X(x)) : = a(x)(X(x)), \/x £ M, 

defines a pairing between Q 1 (M) and A 1 (M). More generally when q in Q q (M) and P in A P (M ) 
are decomposable, so q = ot\ A • • • A a q and P — X i A • • • A X p , for <x t in f2 1 (M) and Xj in 
A 1 (M), we set 


(q, P) := («i A • • • A a q , X 1 A • • • A X p ) 


0 if P + q, 

det((«i, Xj)) if p = q- 


Since the value of {q, P) at a point only depends on the value of rj and P at this point, we can 
extend by bilinearity, in a unique way, this pairing to arbitrary elements q in Q(A/j and P in 
A(M). Furthermore, it is easy to check that if q is in Q P (M), then 

{q, X 1 A • • • A X p ) = q(X u ...,X p ). 

If X is a vector field on M, the inner product i(X) is a derivation of degree —1 on the graded 
algebra Q(M) and since the exterior derivative d is a derivation on A(M) of degree 1, it follows 
that the Lie derivative with respect to A", given by Cartan’s formula: 

Cx ■= [i(AT), d] = i(X) od + do i(X), 

where [•, •] means graded commutator, is a graded derivation of degree 0 on il(M '). C x can also be 
defined on A(M) as the unique derivation of degree 0 such that, for / in A°(M) and Y in A 1 (M), 

C x f = X(f) and £ x Y = [X,Y], 

where [A", Y] is the usual Lie bracket on vector fields. Furthermore the Schouten-Nijenhuis bracket 
is defined as a natural extension of the Lie derivative with respect to a vector field on A(M). More 
specifically, it is defined as the unique bilinear map [•, -]sn : A(M) x A(M) —> A(M) such that, 
for / and g in A°(M) = X(M), X £ A 1 (M ), P £ A P (M),Q £ A q (M) and R £ A r (M), 

a_ ) [f,g]sN = 0 

b-) [A", Q]sn = £xQ 

C-) [P,Q]sn = -{-l)^ q - 1) [Q,P] S N 

d-) [P, Q A R] sn = [P, Q]sn a R+ (-1 )^ q Q A [P, R\ S n 
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From these properties it readily follows that [ P , Q\sn belongs to A p+q 1 (M), therefore the last 
property means that the endomorphism d P : A(M) —> A(M) given by 

dpQ := [P,Q\sn, (1) 

is a derivation of A(M ) of degree p — 1. A somewhat long, but otherwise easy, induction, based 
on the defining properties, gives 

(_!)■>-!)(—1) [p [Q, p] SK ] SJV + (_1)<«-D0.-1) [Q, [p P) s p SM 

+ [p,qM sw = o, (2) 


which is called graded Jacobi identity. 

This, together with bilinearity, c-), and the fact that [ P , Q\sn belongs to A P+9 ~ 1 (M) means that 
A(M), equipped with the Schouten-Nijenhuis bracket, is a graded Lie algebra when the degree of 
P in A P (M ) is declared to be p — 1, not p. So, for instance, vector fields would be the homoge¬ 
neous elements of degree 0 under this new grading of A(M). To perform computations with the 
Schouten-Nijenhuis bracket it is convenient to extend the definition of the interior product. If r/ is 
in f2(M), / is a function and X 1 ,... ,X p are vector fields we set 

i(f)v ■= f d and i(Xi A • • • A X p )r) i(Xx) o • • • o i(X p )p. 

In general, i(P) is defined in such a way that the map i(-): A(M) —> £ , where £ (f2(M)) 

is the space of endomorphism of £l(M), is P(M)- linear. In particular, i(P A Q)rj = i(P){i(Q)rj), 
for all P and Q in A(M). Furthermore, when q is a p-form 

i(X r A • • • A X p )q = i(X,) o ■ ■ ■ o i(X p )rj = V (X p ,..., X,) = (-1 )^rj(X 1 , ...,X p % 
therefore for any P e A P (M) 

i{P)r ] =(-l)^ E {r ] ,P). (3) 

Unfortunately i(P), in general, is not a derivation, which complicates computations. Nevertheless, 
another simple induction gives 


K[P,Q]sn) = [[i(P),d\,i(Q)], (4) 

where the brackets on the right are the usual brackets on the algebra of endomorphisms of A(M). 
Notice that when P = X is a vector field this reduces to the well-known relation among interior 
products and Lie derivatives: i(CxQ) = [Cx,i{Q)\- 
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If P is ap-vector and 77 is a (p — l)-form, then i(P)rj = 0 and i(P) o i(f)rj = i(P)(fr /) = 0. 
These, together with © and © entail 

(V, [p, f} SN ) = (-1 f] SN )r, 

= (-i )*^[[i{P)AAQ)]v 

— (— 1)~ ^ “ (i(P) ° d o i(f)r) — (—1 ) p d o i(P) o i(f)r) 

— i(f) o i[P) o dr] — (— 1 ) p i(f) ° do i(P)r^J 
= (-1) (P 2) 2 (p 1} (i(P) o do i(f)rj - i(f) O i(P) O dry) 

(p—2)(p—1) 

— (—1) 2 i{P)(df Ary) 

= (_l)( P -i) 2 (d/A?/; P) 

= (—l)^ -1 )^ -2 ) (ry A df, P) 

= (r i Adf,P). 


Repited use of this gives 

(df i A • • ■ A df p , P) = (dfi A • • • A d/ p _ l5 [P, / p ] S jv) 

■ • [[-P, /p]sW, /p-l] SJ v’ 


,/l 


SJV 


(5) 


3 Poisson Manifolds 


A Poisson structure on M is a skew-symmetric R-bilinear map {■,•}: P(M) x J~(M) —* P(M) 
satisfying the Jacobi identity: 


{/, {d, Ml + 1A {/, d}} + {d, {h, /}} = 0, V/, 0, h G P(M), 


and such that the map Xf = {•, /} is a derivation of the associative and commutative algebra 
T(M), for each / G P(M), or in other words, Xf is a vector field, usually called a Hamiltonian 
vector field, and / is said to be the Hamiltonian of Xf. This property characterizing derivations of 
the associative and commutative algebra P(M), {<7i<y 2 , /} = < 7 i{< 72 , f}+ 92 {di, /}, called Leibniz’ 
rule, is very important and gives a compatibility condition of the associative and commutative 
algebra structure in P(M) with the Lie algebra given in P(M) by the Poisson bracket. 

To construct Poisson structures let A be an element of A 2 (M), if / G A°(M) and g G A°(M) 
are two functions, using ©, we define a third function by 

{/, 9} ■= A (df, dg) = -A (dg, df) = -(dg A df, A) = - [[A, f) SN , g] SN • (6) 

By construction Xf := [A, f]sN is a vector field, and the defining property b-) entails 


Xf(g) — Cx { g — [Xf,g]sN — ~{fig} — {9i /}• 

In particular 

{9i{hif}}= [A, g]sN, [[A, h]sN,f] SN 


(V) 


J SN 


- [ X gi i X hi /]stv] SJV - £x g o Lx h f- 
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By the same token {h, {/, g}} = —{h, {g. /}} = —C Xh ° £-x g f • On the other hand, the graded 
Jacobi identity ©, the defining property b-), and © entail 


{/, {9 , h}} = -{{ 5 , h}, /} = {[X g , h} SN , /} = - [A, [X g , h] SN ] s , f 


SN 


[Xg, [h, A] SA t] 57V , / 


SN 


= -[[X g ,X h ] SN ,f] 


SN 


[h, [A ,X g ] SN ] SN ,f 
[h, [A,X 5 ] 5jv ] sjv , / 


SN 


SN 


£'[Xg,X h \f \]^i IA, Agjgjv] f 


SN 


Altogether gives 


{/, { 9 , h}} + {g, {h, /}} + {h, {/, g}} 


Furthermore, from the graded Jacobi identity 


(• £x g o C Xh - Cx h ° Cx g - C[ X : g ,x h ]jf 
~ \h,[A,Xg\ SN ] SN ,f. sN 
- [h, [A, X g ]sN] SN , f 1 ■ 


0 = [A, [A, g] 5 Ar] SN + [ A, [A, g]sN] SN + [ 9 , [A, A \ SN ] SN 
= 2 [A, [A, g] SN ] SN + [g, [A, X\ SN ] sn 
= 2[A, X g ]sN + [g, [A, A] SA r] SJV , 

therefore, by © 


{/, { 9 , h}} + {g, {h, /}} + {h, {/, g}} 


h, [g, [A, A] 5 Ar] SAr 

CO 

1 - 1 

h, [[A, A] 5 Ar,^] SJV 

,/ 

SN 


J SN 


SN 


[[A, A]stv, g] SN : h 


SN 


J 


SN 


-~{df A dh A dg , [A, A \ SN ). 


( 8 ) 


(9) 


It follows that the bracket defined via a bivector field A is a Poisson structure exactly when 
[A,A] SA r = 0, and when this happens we say that A is a Poisson tensor. This elementary, but 
clever, computation was first performed by Lichnerowicz in lF70l . who also realized that, when A 
is a Poisson tensor and P is in A P (M), the graded Jacobi identity implies 


0 - (-l) p P [A, [A,-P]sw] SA r - [A, [p, A]sw] SJV + (-l) p ' [p, [A, A] sw ] SAr 
= 2(-ir 1 [A,[A,F] OT ] SK . 

In other words, for Poisson tensors, the derivation d A , defined as in (Q]) by d A P = [A, P] SN, 
satisfies the cocycle condition 

d\ 0 d\ = 0. 
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On the other hand, from ([8]) we see that, for Poisson tensors, [Xf. A]sx = 0, this together with the 
graded Jacobi identity and © give 


[ x fi x g] ~ [ x f> - ~[A, [g,X f ] SN ] SN + [g, [X f , A] 5 iv] SiV 

= [A, l x fi9]sN] SN = -[A, {f,g}] 

= *-{/*}■ do) 


The converse is also true: a Poisson structure determines a Poisson tensor. To see this let £ a 
denote a set of local coordinates on M, then, using the summation index convention, 


x f = x f(Q 


d_ 

di< 


{£a,/} 



hence 


{f,g} = X 9 (f) 


{£a,g} 


df_ 

d£ a ' 


Thus, 

dq 

= -{g,£a} = -{6>, £a}^r- = {£a, 6>} 

and the local coordinate expression of the Poisson Bracket becomes 


dg 

Wb’ 


{f,g} = {&»,&} 


dg df 
d£b d£ a ' 


(ID 


( 12 ) 


Therefore to compute the Poisson bracket of any pair of functions is enough to know the funda¬ 
mental Poisson brackets 

Aaft £&}• 

Moreover, the value of {f,g} at a point m £ M does not depend neither on / nor on g but 
only on df and dg, as explicitly shown in (fl2l) . hence from the Poisson structure we get a twice 
contravariant skew-symmetric tensor 


A (df,dg) := {f,g}. 


Indeed, the local coordinate expression of A is 


A 


a d 
A ab A 

dfa 


d_ 

Wb 


and if £ = o(f ) is another set of local coordinates on M, then, 


A ab ~ 


{£(1) f,b\ {(fai fb\ {£c; id} 


dfg d(j) b 
dfc did 


A-cd 


d(pa dfb 
die di d ’ 


so the components of A do change like the local coordinates of a skew-symmetric twice con¬ 
travariant tensor which, by I©, it is a Poisson tensor. We are using the convention that in the local 
expression of the wedge product only summands whose subindex on the left hand side term is 
smaller than the subindex on the right hand side term do appear. 
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For any function h G J-’(M) the integral curves of the dynamical vector field X h are precisely 
determined by the solutions of the system of differential equations 


d£a 

dt 


{£a, h} , 


and the dynamical evolution of a function / on M is given by 


dj_ 

dt 


{f,h }, 


(13) 


or in local coordinates 

(]£ _ . df_ dh_ 
dt ab df a df b 

Interesting examples of Poisson manifolds are symplectic manifolds. A symplectic form u on 
M determines a bundle map : TM —> T*M over the identity, which gives rise to the corre¬ 
sponding linear map between their spaces of sections, defined by 


(oAa'))(y) :=(o/pf),y)= w (x,y). 


Since u is non-degenerate c J° is actually an isomorphism; we denote X the inverse map and define 
a bivector A by 

A(a,/3) := w(a; a (a),a; # (/3)), 

if a and 3 are 1-forms. When a = df is exact the corresponding vector field is denoted by 
Xj := X(df), and we say Xj is the vector field associated to f with respect to u. It is actually 
defined by the equation i(Xf)u> = df. Furthermore, let {•, •} be the bracket associated to A via ©. 
Then the vector field associated to / with respect to c o is also the Hamiltonian vector field given 
in ©, explaining why we use the same notation. If £ a , a = 1 ,... ,m, denote local coordinates 
and d/d^i, ..., d/df m , and df \,.... df rn are, respectively, the local basis of A 1 (M ) and O' (M) 
associated to £ a , let B = (B ab ) be the matrix of the linear map relative to these bases, and 
lu = uj a h df,, A dfh the local expression of the symplectic form, then 



Thus, B = (c j a b), and the matrix of associated to these bases is the inverse of B. On the other 
hand, 


du>{X f ,X g ,X h )=Xf(u>{X g ,X h ))+X g (u>{X h ,X f ))+X h (u>{X f ,X g )) 

- u ([X f , X g ],X h ) - w([X g , X h ],X f ) - u([X h , X f ],X g ). 

Nevertheless, u(X g , X h ) = A (dg, dh ) = {g, h }, so Q entails 

X,(u(X„X K )) = X,({g,h}) = ~{f,{g,h }}. 

Also, by 03, [X g , ,Y,,] = [X g , v = .Y ^ /h] , hence 

u([Xg, Xh], X,) = u{X-{ g ,h), X,) = -{{<?, h },/} = {/, {g, A}}- 
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It follows that 


0 = du>(X f , X g , X h ) = -2({/, {g, h}} + {g, {h, /}} + {h, {/, g}}), (14) 

so A is a Poisson tensor. 

Reciprocally, from a Poisson tensor A we get a bundle map A 11 : T*M —> TM, defined by 

(a,A»(/?)) A(a, 0 ). 

In general A* is not a bundle isomorphism. We say the Poisson structure is regular when that is the 
case, and then we denote the inverse map by A b . Thus we have an identification of T X M and T*M , 
for each point x of M and, therefore, an identification of higher order contravariant and covariant 
tensors. In particular, we define a 2-form c o by 

u(X,Y) :=A(A l (A').A l (y)). 

Notice that, from this point of view, Xf = A \df). Indeed, from 0 

(dg, X,) = X,(g) = {,g,f } = A (dg,df) = (,ig,\\df)). 


Therefore, 

u(X f ,X g ) = A(A b (A»(rf/)), A l (A«(<ig))) = A (d/. dg). 

In particular the brackets associated to oj and A coincide. Since the Poisson bracket satisfies the 
Jacobi identity, (fl4l) entails, doj(Xj, X g , X h ) = 0. Given that locally one can consider a ba¬ 
sis consisting of Hamiltonian vector fields, we conclude that to is a closed 2-form, moreover by 
definition it is non-degenerate, hence u, so defined, is indeed a symplectic form. Using local co¬ 
ordinates as above, we see that the matrix of A 11 with respect to the standard bases is (A a6 ), where 
A = A ab ^|- A ^ is the local expression of A. 

4 Souriau’s prescription 

As far as we know, Souriau ll37l l38l l39l was the first to realize that since a dynamical system has 
two pieces, the symplectic form and the Hamiltonian, the equations of motion can be described 
by different data, modifying one or the other component. Thus, a perturbed dynamical system can 
be described starting from the free case by modifying the Hamiltonian, as was classically done, 
or simply by changing the symplectic form (see also fTTlO . This idea of adding an extra term to 
the symplectic form was successfully exploited by Souriau in his study of the orbit method, and it 
is what is behind the exotic mechanics, and several other models where non-commutativity of the 
variables is employed. But before we tackle that, let us consider a more down to earth example. 

The classical method to derive the Lorentz equations in a relativistically invariant way is to use 
the so called minimal coupling , which consists in making the substitution p ha p — eA inside the 
free Hamiltonian, where A is the vector potential of the electromagnetic field and e is the electric 
charge. Thus, the starting point is the cotangent bundle T*M, of a manifold M, endowed with its 
canonical symplectic form cu 0 = d0 0 , where 8 0 is the canonical 1-form given, in local cotangent 
bundle coordinates (qi,Pi), induced from local coordinates (q,) on M, by 8 0 = p t d/p, together with 
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a Hamiltonian H : T*M — » M, which one replaces by H A = H o 0” 1 , where (j) A ■ T*M —>■ T*M 
is the bundle map over the identity given by cj) A (q,p) = {q,p + eA(g)), and A = A(<?) dq t is 
a basic 1-form on T*M. The Hamiltonian vector field Xh a associated to this new Hamiltonian 
H a = that leads to the equation of motion, is given by 

i{X HA )de 0 = dH A = {<j>2 1 Y{dH). 

On the other hand, and with an abuse of notation, we denote cr both the 1-form on M defined by 
a = eA l (q) dqi, as well as the basic 1-form on T*M obtained by pulling back a by the canonical 
projection n : T*M —> M. Then, 


(j)* A (d0 o ) = (p* A (dpi A dqi) = d4>* A (Pi ) A d(fi* A (qi) = d(pt + eAi(q))Adqi = dd 0 + da, 


and since 0 is a diffeomorphism, 

i((j>A*(X HA ))(d9 o + da) = i(<f>A,{X HA ))(<j>* A d6 o) = <j>* A [i{X HA )de Q ) = 0a((0a) _1 ( dH )) = dH. 

In other words, by adding the extra term da to the symplectic form, which, by the way, it is a basic 
2-form (i.e. locally it is of the form g tJ (q) dqi A dqj, so it does not involve the p’s), we see that 
the vector field 4>a*(X H a ) is the Hamiltonian vector field associated to the original Hamiltonian H 
with respect to this new symplectic form, and we obtain the same equations of motion. 

If u = + -ppj dqi A dqj, where gij(q,p) is a skew-symmetric matrix, and the Hamiltonian 

vector field X H is X H = V,d (H + W,d Pi , then 

i{X H )u = (■ i(X H )dqi ) A dp t - dq t A (■ i(X H )dpi ) + ^gij(i{X H )dq^) A c% 

- ^9ij dqi A (i(X H )dqj) 

= ~Wi dqi + Vi dpi + g t j Vj dqi 


The equation i(X H )u> = dH, entails 


T r dH J TI7 OH dH 
Vi = -a— and W t = , 

OPi OQi OPk 

therefore by (|71) the Poisson bracket associated to u is given by 

{F, H} = J- - + gM- = {F, H} o + 9i ™L = X H (F), 


dpi dqi dqi dpi dp { dpj 


dpi dpj 


where {-,-} 0 stands for the Poisson bracket corresponding to u 0 . It follows that the generalized 
(Hamiltonian) vector field is 

_ dH d 

— A H + gija ;—• 

opi dpi 


The equations of motion are given by 











Our construction can be extended easily to a more general Souriau form 


u = u 0 + \gij dqi A dqj + \fij dpi A <%, 


and the equations of motion are then given by 

dq k _ Oil | f Oil dp k dH 8H 

7. r\ Jki r\ 7 j i r\ ” 1 ” Qik r\ 

dt opk dt OQk dpi 

Then if we assume the Hamiltonian is of the form 


H (q,p) = ^^ + V{q), 

with the potential energy V depending only on the configuration coordinates g*, the equations of 
motion are 

dq k _Pk f 9V dp k _dV Pi 

dt m oqi dt oq k m 

These are equivalent to the modified Newton’s second law £30! 4.5, 501I5TI 


m 


d 2 q k 
dt 2 


dV 

dqk 


pi d 

+ 9ik -h m— 

m dt 



(15) 


The second term of equation (fl5l) is a correction due to the noncommutativity of momenta and the 
third term is that of noncommutativity of coordinates. It follows that even for the case V = 0 the 
particle accelerates because of the noncommutativity of momenta. 

The second procedure has the advantage that it works even when only the 2-form is globally 
defined, with no reference to the 1-form 9 0 made. In lT72ll this idea was generalized to the case 
of a classical particle in the presence of a Yang-Mills field. When the Poisson manifold M is 
the tangent bundle TQ of a //-dimensional manifold, so M is known as the velocity phase space , 
Souriau also proposed to describe the dynamics not on phase space but in what he called evolution 
space , with coordinates (xi, Xj,t). His idea was to join the symplectic form u on phase space 
with the Hamiltonian by considering the two-form to — dH A dt on the evolution space, and then 
perform the minimal coupling recipe. This allows to recover the Euler-Lagrange equations, and 
it is equivalent to Faddeev-Jackiw construction If29ll731 . Recently, Bolsinov and Jovanovic lf74l 
considered G'-invariant magnetic geodesic flows on coadjoint orbits of a compact Lie group G, 
where o is the Kirillov-Kostant two-form. 


4.1 Souriau’s formalism and exotic mechanics 

For concreteness let us ponder the ‘exotic’ plane studied by Horvathy in ETl . Thus, we consider 
the dynamical system (T*M 2 , H 0 ), where iieR, 


c 0 $ = dqi A dpi + dq 2 A dp 2 — id dpi A dp 2 


and 


H 0 = 


pl + pj 

2m 
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The 2-form u# is not only closed but exact, and as the associated map iJ# : A 1 ( T*R 2 ) —» fl 1 (T*M 2 ) 
is given by the matrix 

/0 0 -1 0 \ 

0 0 0 -1 
10 0 '!? 

\0 1 -!? 0/ 

which is regular for any value of d, the 2-form u)$ is symplectic. The inverse matrix is 


u# = 


= 


( 0 

d 

1 

°\ 

-1? 

0 

0 

1 

-1 

0 

0 

0 

l 0 

-1 

0 

0/ 


which corresponds to the Poisson structure associated to the bi-vector 

. d d d d d d 

-A-i? = d—— A —-h -— A -h -— A -— , 

oqi oq 2 oqi dpi uq 2 op 2 


so the fundamental Poisson commutators are 


{91,92}=$, { 9 i, Pi} = {92, P2} = 1 , {91, P2} = {92, Pi} = 0 , (pi,p 2 } = 0 , 

and the dynamical vector field is given by 

9 i = ( 9 i, #0} = —, 92 = {92, #0} = —, Pi = (Pi, H 0 } = 0 , p 2 = (p 2 , H o } = 0 . 

m m 

In this way we obtain a 1-parameter family of symplectic structures for the free particle. 

These symplectic structures are the sum of two symplectic structures that are invariant under 
rotations in the plane, one of them being the canonical symplectic structure on T*M 2 . The gener¬ 
ating function of such 1-parameter group will be the function / satisfying 

i{X)u# = df, 


where X is the vector field that is the cotangent lift of the rotation generator in configuration space, 

d d d d 

X = 9i~-92 "7 T Pi “7-P 2 7 ^ • 

oq 2 oq\ OP2 Opi 

Since 

i{X)u# = p 2 dqi - pi dq 2 - 92 dpi + 9 i dp 2 + $(pi dpi + p 2 dp 2 ), 

we find that the generating function is given by 

d 


/( 9 i, 92 ,Pi,P 2 ) = 91P2 -92P1 + 2 (p? +P2) • 


We now apply Souriau’s minimal coupling procedure, so we introduce a basic 2-form a = 
B(qi, q 2 ) dqi A dq 2 , which is closed and can be interpreted as a magnetic field, and consider the 
closed 2-form 

6c?$ 5<j .— 7T O' • 
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The corresponding linear map uj^ a : A 1 (F*Wl 2 ) —> Q 1 (T*M 2 ) is represented by the matrix 

/ 0 5 -1 0\ 

b _ -5 0 0-1 

^ - 10 0 0 ’ 

\ 0 1 -0 0 / 

whose determinant is (1 — 0 5) 2 , therefore u^ t(T is regular in the points where BO ^ 1; the inverse 
matrix being given by 


( 0 

0 

1 

°\ 

-0 

0 

0 

1 

-1 

0 

0 

5 

l 0 

-1 

-5 

0 / 


which corresponds to the Poisson structure defined by the bi-vector 

l/ <9 d d d d d d d \ 
’ ?,cr 1 — 0 B \ dqi ^ dq 2 dqi ^ dp\ dq 2 ^ dp 2 dp\ ^ dp 2 ) ' 


The corresponding fundamental Poisson brackets are 


{gi,®} = 

{Pl,P2> = 


0 


1-05’ 

B 

1 -OB 1 


{9l, Pi} = { 02 ,P 2 } = 


1-05’ 
{91,52} = {92, Pi} = 0. 


When the Hamiltonian is 


H = 7T- + v (9i, 92), 
2 m 


using (HU), (H3b and (fl6l) . we see that the time evolution is given by 

0 dV pi . 1 dV 


91 = 

92 = 


1-05 <9g 2 m(l ~‘0B) ’ 

0 <9H p 2 


1 — 05 dqi m( 1 — OB) 1 


Pi = 

P2 = 


+ 


5p 2 


1 — 0 B dq 2 m(l -05) 


(16) 


1-05 dq x m( 1 — 05)’ 
1 <9H Bpi 


The system is still invariant under rotations if 5 is a rotationally invariant function, i.e. 5 is a 
function of qf + q 2 , B = B(qj + q|), and the generating function for the infinitesimal generator of 
rotations is 

$ 1 

/( 9 i, 92 ,Pi,P 2 ) = 91P2 -92P1 + 2(01+02) + ~B(ql + q%). 

On the other hand, when B is constant, and 50 = 1, the determinant of c J°g a is zero and the rank 
of uj\ a is two, the kernel of the 2-form uj$ jCT being generated by the vector fields 


Xi 


d d 

0-h - 

dq 2 dpi ’ 


d d 

and X 2 = -0--h — . 

dqi dp 2 


The solutions of X\F = X 2 F = 0 are to be found from the method of characteristics and turn out 
to be the functions which depend on £1 = qi + 0 p 2 and £ 2 = q 2 — 0 p x . This suggests the change of 
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variables (gy, g 2 ,Pi,P 2 ) ^ (6, &,Pi,P 2 ), i.e. <71 = Q -$£> 2 , <?2 = 6 + ^Pi- In such coordinates, 
A"i = d/dpi, X 2 = d/dp 2 and becomes 

cc^ ]Cr = —5 i A • 

This means that the quotient manifold T*Wi 2 / kertu^ i(T is parametrized by Q and £ 2 which, more¬ 
over, are Darboux coordinates for such 2-dimensional symplectic manifold. 

Using the same idea with commutators Nair and Polychronakos If32l described quantum me¬ 
chanics for both the non-commutative plane and the non-commutative sphere, and proved that the 
Landau problem for the non-commutative plane can be recovered as the limit of large radius of the 
Landau problem for the non-commutative sphere. 


4.2 Nonrelativistic anyon model in Souriau formalism 

Let (qi,q 2 ) be orthogonal Cartesian coordinates in Q — 1 
IF{TQ) of the free particle 


\ and consider the Lagrangian L 0 in 


L'oiqi, < 72 , Vi,v 2 ) = -m{v\ + vl) . 

Let M be the graph of the corresponding Legendre transformation. This is the submanifold of the 
Pontryagin bundle TQ © T*Q given by the constraint functions 

\ i (q 1 ,q 2 ,v 1 ,v 2 ,p 1 ,p 2 ) — Pi — mvi . 

Let fiGlbea constant and consider TQ endowed with the exact 2-form ui defined by 

Ui(qi,q2,v 1 ,v 2 ) : = ndvi A dv 2 . 

In the spirit of Souriau’s idea, we consider the closed 2-form on TQ © T*Q 

u : = —pr*aq + pr 2 cc 0 , 

where pr x andpr 2 are the natural projections pr x : TQ®T*Q —> TQ andpr 2 : TQ®T*Q —> T*Q, 
tc 0 is the canonical symplectic structure in T*() and oj\ e f 1 2 (TQ) is as before. The corresponding 
map c o° : A 2 (TQ © T*Q ) — > VI 1 (TQ © T*Q ) is represented by the matrix 


c o = 


/0 0 
0 0 
0 0 
0 0 
1 0 
\0 1 


0 0-1 

0 0 

0 K 

-K 0 0 

0 0 0 

0 0 0 


0\ 
0 -1 
0 0 
0 
0 

0 / 


which is regular for k ^ 0. In this case, the 2-form oj is symplectic, and since the inverse matrix is 

/ 0 0 0 0 1 0\ 

0 0 0 0 0 1 

0 0 0 -1/k 0 0 

0 0 1/k 0 0 0 

-1 0 0 0 0 0 

\ 0 -1 0 0 0 0 / 


oj* = 
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we obtain the bi-vector field on TQ © T*Q 

A a ^ <9 a 1 c> a ^ 

dqi dpi dq 2 dp 2 n dv± dv 2 ’ 

and the corresponding Poisson structure is defined by the following fundamental relations: 

{Qi,Pi} = U2,P2} = 1 , {91,92} = {pi,M = {91, M = {Pi, 92} = 0, 

{u 1 ,u 2 } = — , {vi,P 2 > = {w 2 ,Pi> = 0 . 

K 

The two constraint functions Ai and A 2 are second class constraints, because 
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Tfi 

{Ai, A 2 } = {pi - mvi,p 2 - mv 2 } = m 2 {v l: v 2 } = -. 

K 

They define a four dimensional symplectic manifold. 


5 Feynman-Dyson’s method and non-commutativity 


In this section we review the Feynman’s derivation of Maxwell’s equations {6), in the framework of 
a tangent bundle, so the Poisson manifold M is the tangent bundle TQ of a configuration space Q. 
In terms of local tangent bundle coordinates in TQ induced from local coordinates in Q, denoted 
Xi and ±i, a general Poisson bracket on TQ is locally given by 


. , dg df dg df dg df r . . , dg df 

{/, 9} = K, + R, x b } — — + {x a , x b } — — + R, x b } — — . (17) 


r\ r\ 1 J r\ • r\ 1 J r\ r\ • 

OX 5 OX a OXOX a OX ft ox a 

The asumptions in [j6] are Newton’s equations of motion 


dx b dx a 


mxj = Fj(x, x), 

i.e. the dynamics is given by the second order differential equation vector field 

d d 

r = Xi ^- h Fi(x,x)—, 

OXi OXi 

together with the fundamental brackets 


{xi,Xj} = 0 and m{xi,Xj} = Sij. (18) 

The goal is to determine the other fundamental Poisson brackets, and as {x tJ x 3 } must be skew 
symmetric it can be written as 

{ X 1 ■ X j } n^ijkFk (x, X) j 

m z 

where £ tJ k denotes the fully skew-symmetric Levi-Civita tensor and B is defined as the magnetic 
field. Now (fl8l) implies that 

t x . _p.t = J 

3 m d±i ’ 
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and using the derivation property for the time derivative of the second equation in (fl8l) . i.e. assum¬ 
ing that the vector field T is a derivation of the Poisson structure, we get 

1 1 

{Xi, Xj} = - {Xi, Fj} = —£ ijk B k (x, x ), (19) 

m m z 

and the Jacobi identity for the functions ay, x 3 . x k . 


{Xi, {Xj, Xfc}} + {x k , {x h Xj}} + {Xj, {x k , Xi}} = 0 


entails 


0 }x r , Bg} 


1 dB s 
m dx r 


In particular B does not depend on the dotted variables. Moreover from the Jacobi identity with 
three different velocities we obtain 

div B = 0, 


which reveals that the flux of the field B through a closed surface is zero, and that magnetic 
monopoles do not exist! 

On the other hand, since B does not depend on x t , equations (fT71) and (fl9l) entail that F is at 
most linear in such variables, therefore we can define another field E, called the electric field, by 
Ej = Fj — £j ki x k Bi. Using repeatedly all the equations above, one arrive to Maxwell’s equation 
corresponding to Faraday’s law of electrodynamics, in the setting suggested at the beginning of 
this section, namely rot E = 0 in the autonomous case, or in general 


<9B 

~dt 


rot E = 0, 


a magnetic field that is changing in time produces a non-conservative electric field. We refer the 
reader to [33, 34j for details. 

To obtain a dynamic different from the standard formalism of electrodynamics, one needs to 
modify the fundamental brackets in (fl8l) . The first idea is to use Souriau’s technique, namely to 
replace the first fundamental bracket by 


{xi,Xj} = gij(x), 


where g i3 is an arbitrary skew-symmetric matrix of functions, fulfilling the constraints that the 
Poisson bracket properties impose, and keep the other assumptions. In particular, the Jacobi iden¬ 
tity 

{Xi, {Xj, x k }} + {x k , {Xi, Xj}} + {Xj, {x k , Xi}} = 0, 


entails 

0 = {x k ,gij} = {xkiXi}^ 1 + = -—^FL. (20) 

oxi oxi m oxk 

Then the matrix g i3 is a constant skew-symmetric 3x3 matrix, which is an interesting, but some¬ 
what restrictive, case. We then modify Souriau’s idea and settle for 


{xi,Xj} = gij(x,x). 
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Accordingly, (l20l) becomes 



This equation clearly relates the part of the Poisson structure on the base (the positions) with the 
part of the Poisson structure on the fibre (the velocities). Hence if the Poisson structure on the base 
is known one can compute the fundamental brackets on the fibre. 

On the other hand, from the Jacobi identity among (x*, Xj,x k ) we obtain 


9jk } T {x kl Qij\ T 9ki\ 0, 


and this leads to another constraint: 



( 21 ) 


Note that if d x denotes the exterior derivative on the vector space T X Q, for x in 0. then the term 
inside the parenthesis in (1211) are the local coordinates of d x u>, where cu x is the 2-form in Q 2 (T X Q) 
defined by u x = g l: j(x, x)dx t A dij. On the other hand, if A is the bivector in A 2 (TQ) given 
by A = (jjj (x, x)-£- A the terms outside the parenthesis in ([Til) are the local coordinates of 

[A, A] sn- Therefore equation (CHT) is fulfilled when 


d x u = 0 and d^A = [A, A]sn = 0, 


that is, when Cj x is closed and A is a Poisson tensor. 

Furthermore, similar ideas as in the commutative case, using the other Jacobi identities, leads, 
see IfDfl . to the modified Gauss law 


div B = —— B • V x B 


m 


Odd 
dxi ’ d±2 ’ 8x3 


where V = ( 


), and also to 



which is what replaces the Maxwell equation corresponding to Faraday’s law. 

5.1 Generalized Lorentz force equations 

Consider the Hamiltonian dynamical system on TM 3 , where the Hamiltonian and the symplectic 
form are given respectively by 



and 



19 











Assume that, in local coordinates, the Hamiltonian vector field is written as X H = Sid Xi + Rjd Xi . 
The equation i(X H )Q = dH becomes 


1 

—Xi = 

m 

1 

—X2 = 
m 


—S i + 921R2 + 931R3, 

m 

— S 2 + <7l2-Rl + .9.32 R 3 , 


m 


—£3 — — S3 + g^Ri + 923-R2, 
m m 

On the other hand, from (fl3l) and ([7]) we obtain 


dxi 

dt 


Si 


and 


dp 

dx\ 

dp 

dx 2 


8 x 3 





+ R '2 S 3 ~ 
+ B 3 S 1 — 
+ B 1 S 2 — 


B 3 S 2 , 
5i5 3 , 
B 2 S!. 


d±i 

dt 


= Ri- 


( 22 ) 


Therefore if we assume that X H is a second order differential equation, namely that S, = x t , then 
= Ri, and the right column of (l 22 b entails 


1 d 2 x 1 

8 <j) 

m dt 2 

dx 1 

1 d 2 x 2 

dtp 

m dt 2 

dx 2 

1 d 2 x 3 

8 <p 


m dt 2 dx 


+ x 3 B 2 - X 2 B 3 , 

+ X1B3 — X3B1 , 
+ x 2 Bi - ±iB 2 , 


whereas the left column provides the constraints 


n ( 

0 = mg 21 -x— 

V ox 2 

n ( 

0 = mg 12 

V ox 1 

n ( d< ^ 

0 = mja \xr l 


X1B3 + x 2 Bi 
X3B2 + X2B3 
X3B2 + X2B3 


+ rrujsi 
+ mg 32 
+ mg 2 3 


(<W>_ 
V<9x 3 
/ 8<t> 

Vc)a:3 

f <90 

\dx 2 


- x 2 Bi + 


- x 2 Bi + 


- X 1 B 3 + 


xiB-^j, 

X3Blj. 


(23) 


In particular when V0 = —eE, (1231) is a generalized Lorentz force: a force experienced by a 
charged particle moving in an electromagnetic field, subject to a system of constraints. 

Another interesting class of systems can be studied via the “generalized Souriau form” 


d ) 0 = d£i A dxi + gijdxi dij. 


It is a mixture of a symplectic and a gradient structure, known as a metriplectic system. The sym¬ 
metric bracket associated to the metric tensor incorporates the dissipative structure of the system. 
The Leibniz vector field X h associated to a function h G C°°(M ) satisfies X h = V/i, i.e. X h 
generates a gradient dynamical system. In local coordinates the vector field Xh is given by 


Xh flij 


dh 8 

dxj dxi 


and the corresponding bracket in this context is called a Leibniz bracket. 
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6 Volume preserving mechanical system related to Souriau form 


Another interesting class of dynamical systems that generalize the Hamiltonian systems, where 
noncommutativity is also possible, was introduced in [|47l 481. 

Let (M, Co'o) be a 2 n-dimensional symplectic manifold, a vector field A" is said to be symplectic 
if fxw 0 = 0, from the Cartan identity, this is equivalent to i(X)loq being closed, in particular 
every Hamiltonial vector field is symplectic. On the other hand, we say a vector field A" preserves 
the volume if Cx^o = 0; here and in what follows powers are meant with respect to the wedge 
product. Since 




s=0 


where <f>/ is the flow of A, it follows that <!>/ does preserves the volume form cjq . Furthermore, a 
simple induction gives 


— kC x ^o A Wq 1 


In particular, we see that every symplectic vector field preserves the volume, but the converse is 
not true in general. The divergence of a vector field X is defined as the unique function div X in 
C°°(M ) such that 


£x^q = div X cjq. 

Therefore X preserve the volume if, and only if, it is divergence free. Let (xi,, x 2n ) be Darboux 


coordinates, then u 0 = dx t A dx n+i , and if A" = Yli=i VA, it is easy to check that 



We now describe a procedure that produces dynamical systems that preserves the volume. 
First consider the map F: A 1 (M) —> H 2n_1 (M) given by F(X) := i(X)u}q. Using Darboux 
coordinates, simple combinatorial arguments entail 








(24) 


where as usual the hat means that the term is to be deleted, in particular 


n(n — 1 ) 

i(X)ujQ = (—1) 2 n\i(X)dxi A • • • A dx 2n 



(25) 


therefore F is surjective. Moreover, if X(p) ^ 0 and F(X) = 0, we can locally find a basis 
A"|, X 2 n, of vector fields such that X 2 = X, then 
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which is absurd since cc 0 is nondegenerate, hence F is injective and therefore a linear isomorphism. 
Since C x u o = di(X)ujQ, under F, the space of volume preserving vector fields corresponds to the 
space of closed (2 n — 1)-forms. Thus, if g is a 2-form 

w-= A^r 2 )) 


is a volume preserving vector field. In particular, if lu is the Souriau form, written in Darboux 
coordinates as 

uj dx n -\-i A d'Xi -(- —(jijXn-\~i A dx n +j T -Bijdxi A dxj , 


then 


d(u; A 0Jq 2 ) = du A Uq 2 


1 dg {j 


2 dxk 
1 dB 


dxk A dx n+ i A dx n+ j A uJq 


n—2 




2 <9x 


n+/c 


dx n+ k A dxi A dxj A cu ( 


n—2 


Now, from 


dxk A dx n+ i A dx n +j A cCq ~ ~ 2 = <5jw dx n+ j A dxk A dx n+ i AcUq 2 

- S kj dx n+i A A A cUq " 2 


Similarly 


Thus 


(-l) n 

n — 1 


(-i)’ 


$ki dx n -\-j A cJq &kj dx n -^-i A 


n— 1 


dx n+ k A dx* A A o;q 2 =-— ( S^j dx^ A o;q 1 — A a/ 1 1 


,n—2\ 


n — 1 

l) n / dg^ 


Jki axj /\ uj 0 


d(u A ay, . . „ 

v 0 ' n - 1 V 5rc. 


dx n -^-j A cJq 


n—1 


dB, 


<9x 


dx, A ca 


71 — 1 


n+/c 


Let X be the vector field 


_ _ dg u d dB u d 

A. C77, ^ ^ 


with 


dx k dxi n dx n+k dx n+ i ’ 

then, using (1251) and (l24l) 

z(X)uJq = nl(— 1) ( 2 )+l_1 c n dx 1 A • • • A dx^ A • • • A dx 2n 


C"n. 


n(n-l) 

Hj 2 

n(n — 1) 


dan 


95, 




<Te. 


n+/c 


dxi A • • • A dx n+ i A ■ ■ ■ A dx : 


2 n 


. n(n — 1 ) . 1 • -1 

— n!(—1) 2 +n+ *“ 1 c 

Qg^ . _ ____ 

— (n — 2) !(—!)” ——- dx n+i A dx i A ■ ■ ■ A dx, A ■ ■ ■ A dx n+i A ■ ■ ■ A dx 2n 


dx 


(n-2)!(-l)' 


dB ki 


dx. 


n+k 


dxi A dx 1 A • • • A dxj A • • • A dx n+i A ■ ■ ■ A dx 2n 


= (n- 2)!(—l) n+( 2 J dx n+i - >? Bk% dx { , 

\ dx k dx n+ k 


A dx 1 A <ir n+1 A • • • A dx, A dx n+i A ■ ■ ■ A dx n A dx 2n 


= d(u A oJo )- 
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In other words, X u = X, so we can associate a volume preserving flow with the Souriau’s form, 
and the equations of motion are given by 

dxj _ (-l) n dg ki 
dt n — 1 dx k 
dxn+j _ _ (~l) n dB ki 

dt 7i — l dx n+k 

A Nambu-Poisson system is a volume preserving flow, on a Nambu-Poisson manifold of order 
2 n, determined by (2 n — 1) Hamiltonian functions Hi, , H 2n -i £ C°°(M) 

dxi 

,...,H2n —l 'Hi) 1 dd i , ... , H 2 n—i, X\ j. 

In general, C Xlll . H2n _J = i H u ■ ■ •, H 2n -i, /} if X Hu ... :H2n _ 1 is a Nambu-Hamiltonian vector 

field. If rj is a Nambu-Poisson tensor, then Takhtajan 11731 proved jC Xi1i h 2 , V — 0. 

7 Conclusion and Outlook 

In this paper we have studied the classical non-commutative mechanical systems using Souriau’s 
method of generalized symplectic form. In particular we have explored a large class of non- 
commutative flows which includes the non-commutative magnetic geodesic flows, non-relativistic 
anyon model |76l 77, 78], generalized Lorentz force equation, etc. Souriau’s formalism allows us 
to study geometrically all these non-commutative dynamical systems in an unified manner. The 
dynamics of these systems boil down to generalized Hamiltonian dynamics where the Poisson 
structure can be complicated functions of phase space coordinates and momenta. However, some 
questions should be addressed in near future. 

At first we must consider the quantization of these classical non-commutative system. There 
is an interesting paper If79l which addresses the connection between non-commutative quantum 
mechanics and Feynman-Dyson’s method. Actually, the generalization of Feynman-Dyson’s idea 
to the quantum world would be an interesting subject to be studied. This would lead to unveil 
the close relation existing between the non-commutative geometry and the geometric phases. 
Quantization of these models could give rise to new physics at some very high energy scale 

m m m (Mi m mi m m m m uni m. 

This method can be extended to other directions also. We can simply generalize this construc¬ 
tion to supersymmetric non-commutative systems Il93l |94|| . In other words, we can try to gener¬ 
alize the Feynman-Dyson’s scheme to supersymmetric framework. This would certainly yield a 
generalization of supersymmetric generalized Hamiltonian dynamics. There is a recent upsurge of 
interest in (2 + 1)-dimensional model [92], [931 with a kind of nonstandard noncommutativity, where 
both coordinates and momenta get deformed commutators. So far most of the papers concern time- 
independent systems, it would be rather challenging to extend this framework to time-dependent 
systems. In a recent note Liang and Jiang [|96l studied the time-dependent harmonic oscillator in 
a background of time-dependent electric and magnetic fields. Recently noncommutative quantum 
mechanics |[79[[97ll98l[99 t T001 is becoming a exciting topic to study, it would be interesting for us 
to investigate this subject using geometrical methods of quantizing noncommutative phase space 
mechanics. 
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Possibly one can study the Helmholtz condition analogous to standard classical dynamics. 
Then the corresponding variational formulation can be used to construct Noether symmetries and 
conserved densities. It is known [144[| that the Helmholtz condition connected to (9-dcformcd Pois¬ 
son system is a third-order time derivative equation. One should analyse carefully all these new 
aspects. 

We can also study the field theoretic Poisson brackets on jet space. This will yield an interesting 
class of partial differential equations. One must try to explore its connection to other branches of 
mechanics and geometry, namely, non-holonomic systems, control theory, Finslerian mechanics, 
Lie algebroid theory, etc. 
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